Adiabaticity Criterion for Moving Vortices in Dilute Bose-Einstein Condensates 
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Considering a moving vortex line in a dilute atomic Bose- 
Einstein condensate within time-dependent Hartree-Fock- 
Bogoliubov-Popov theory, we derive a criterion for the quasi- 
particle excitations to follow the vortex core rigidly. The as- 
sumption of adiabaticity, which is crucial for the validity of 
the stationary self-consistent theories in describing such time- 
dependent phenomena, is shown to imply a stringent criterion 
for the velocity of the vortex line. Furthermore, this condi- 
tion is shown to be violated in the recent vortex precession 
experiments. 
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Since the first experimental realizations of Bose- 
Einstein condensation in dilute, harmonically trapped 
atomic gases pL there has been great interest to inves- 
tigate the superfluid properties of these unique quantum 
fluids. Due to the inherent connections between quan- 
tized vorticity and superfluidity, this interest culminated 
as the creation of vortices in trapped condensates was 
demonstrated [Q]. The recent experimental advances in 
manipulating vortices and observing their dynamics are 
providing efficient tools to study the physics of these in- 
teracting many-particle systems and to relate it to the 
quantitative predictions of thermal field theories. 

The structure and, in particular, the stability of vor- 
tices in dilute Bose-Einstein condensates (BECs) has 
been under an extensive theoretical analysis ||. The 
majority of the studies has been carried out within 
the zero-temperature mean-field formalism consisting of 
the Gross-Pitaevskii (GP) and Bogoliubov equations. 
Within the Bogoliubov approximation (BA), the excita- 
tion spectra of vortex states in statically trapped conden- 
sates have been shown to contain at least one mode with 
positive norm but negative energy Q. These anoma- 
lous modes have crucial consequences for the superfluid 
properties of the condensates, since they imply the vor- 
tices to be energetically unstable in nonrotating traps. 
Furthermore, these states have been shown to manifest 
themselves in the precession of the vortex line about the 
symmetry axis of the trap, with the precession frequency 
and direction determined by the excitation energy — 
especially, the negative energies imply precession in the 
direction of the condensate flow |5j . 

The predictions of the Bogoliubov approximation agree 
well with the experiments. The critical trap rotation fre- 
quencies for vortex nucleation can be understood theo- 



retically to good accuracy Q. Also, the precession of vor- 
tices predicted by the GP equation has been experimen- 
tally observed [Q . The precession frequency and, in par- 
ticular, its direction are in line with expectations based 
on the BA. In general, the mean-field theory has turned 
out to be remarkably successful in describing trapped 
BECs, including the vortex states and their dynamics 

I- 

However, the situation changes when the analy- 
sis is taken beyond the zero-temperature BA by self- 
consistently including the effects of the thermal gas com- 
ponent. Stationary self-consistent solutions for vortex 
states within the Popov approximation (PA) and its re- 
cently proposed extensions contain no anomalous modes 
even in the zero-temperature limit [PUlOfl- This is due 
to the partial filling of the vortex core with the noncon- 
densate, which serves to lift the anomalous quasiparticle 
states to positive energies. The positive precession mode 
energies, in turn, imply vortex precession opposite to the 
condensate flow, in evident contradiction with the ex- 
perimental observations and the predictions of the BA. 
In the light of the success of the self-consistent approxi- 
mations in predicting excitation spectra for irrotational 
condensates |ff|,^2[ , this discrepancy is surprising. In ad- 
dition, the close agreement of the BA with the results of 
the vortex precession experiments implies that the mean- 
field approximation itself is not the cause of the failure 
of the PA. 

We suggest that the apparent disagreement between 
the PA and the experiments could be due to incomplete 
thermalization and/or inadequacy of the quasi-static for- 
malism in describing moving vortices. In order to clarify 
the latter possibility, we show in this paper that the va- 
lidity of the quasi-static self-consistent mean-field treat- 
ment in modelling moving vortices imposes for the vortex 
velocity a stringent criterion, which seems to be violated 
in the precession observations reported so far. This im- 
plies that at the observed velocities the quasiparticles can 
not follow the vortex core rigidly, and its structure and 
spectrum are deformed from those of a static vortex. 

In order to describe the dynamics of trapped BECs, 
we use a time-dependent mean-field formalism based on 
the Popov approximation [ p^[ . Working in the grand- 
canonical formalism, we start with the Heisenberg equa- 
tion of motion 



_ d , . . 



H (r)i>(r, t) + gtf(r, t)^(r, t)V(r, t) (I) 
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for the field operator ip(r, t) of a dilute boson gas. Above, 
7Yo(r) = — h 2 V 2 /2m + V tI (r) — fi is the grand-canonical 
one-particle Hamiltonian corresponding to the trapping 
potential 14 r (r) and the chemical potential and the 
coupling constant g is related to the s-wave scattering 
length a by g = AirTi 2 a/m. Inserting into the nonequilib- 
rium average of Eq. (Q) the Bogoliubov decomposition 



ip(r,t) = $(r,t) +4>(r,t) 



(2) 



of the field operator in terms of the c-number condensate 
wavefunction $(r, t) = (ip(r,t)) and the noncondensate 
field operator tp(r,t), and treating the expectation val- 
ues of noncondensate operator products according to the 
Popov mean-field scheme, we arrive at the generalized 
GP equation 

»&^*(r, t) = £(r, t)$(r, t) - gn c {r, t)$(r, t) (3) 

for the condensate wavefunction. Above, C(r, t) = 
Ho(r) + 2gn(r,t), and 



n c (r,t) = |$(r,i)| 2 , 
h(r,t) = (i>\v,t)i>(r,t)}, 
n(r, £) = n c (r, f) + n(r, t) 



(4a) 
(4b) 
(4c) 



denote the condensate, noncondensate, and total den- 
sities, respectively. Correspondingly, within the Popov 
mean-field approximation, one finds 

ih^{r, t) = C(r, t)$(r, t) + g<S> 2 (r, t)ft{r, t) (5) 

for the equation of motion of the noncondensate field op- 
erator. Substituting into Eq. (|^) the Bogoliubov trans- 
formation 



(6) 



of the field operator in terms of the bosonic quasiparticle 
operators a n and cv n , we find that the quasiparticle am- 
plitudes it„(r, t) and v n (r,t) satisfy the time-dependent 
Hartree-Fock-Bogoliubov-Popov (TDHFBP) equations 

d 

ih—u n (r, t) = £(r, t)u n {r, t) - 5$ 2 (r, t)v n (r, t) (7a) 
v n (r,t) = L(r,t)v n (r,t) - g$* 2 (r,t)u n (r,t). (7b) 



Introducing the matrix notations 

' u n (v,t) 



/n(r,*) 



v n (r,t) J ' 



C(r,t) 



£(r,t) -g<£ 2 (r,t) 
g^>* 2 (r,t) -C(r,t) 



(8) 



(9) 



the quasiparticle equations can be expressed in the com- 
pact form 



ift^/„(r,t) = 0(r,f)/„(r,t). 



(10) 



For convenience, we also define positive- and negative- 
sign scalar products of quasiparticle states by setting 



</<l/i> = 



dr K(r) Uj (r) ± v*(r)vj(r)]; 



(11) 



here and henceforth, we suppress the arguments of func- 
tions when they are not needed for clarity. The re- 
quirement that the quasiparticle operators a n , sat- 
isfy canonical bosonic commutation relations implies for 
the quasiparticle states the normalization (fi\fj)_ = 
correspondingly, only states with positive norm are to 
be included in the Bogoliubov transformation of Eq. (|^). 
This normalization can be straightforwardly verified to 
be consistent with the TDHFBP equations. 

In case the mean fields and, hence, the operator C(r, t) 
vary slowly in time, we expect that the solutions of the 
TDHFBP equations may be approximated by solving at 
each instant of time the corresponding quasi-stationary 
eigenequations 



(12) 



This adiabatic approximation is accurate if the transi- 
tion rates, as determined by the exact time development, 
of the quasi-stationary states to each other are negligi- 
ble. In order to formulate this criterion quantitatively, 
we follow the treatment of Ref. JO]. Let {/f°' > } be a 
complete set of solutions of Eq. (|12D; especially, it con- 
tains the zero-energy solution fjf^ oc ($0j ®o) T > where 
$o is the solution of the stationary GP equation. We 
orthonormalize the solutions by requiring 

K/J 0) |/j 0) >_l = *« (^0); </ri/r>+-l- (13a) 
In addition, we may impose the condition |15| 



(/ri// o) ) + =o (<^o). 



(13b) 



In order to analyse the transitions of the quasi-stationary 
states to each other, we expand the solutions of the TD- 
HFBP equation in terms of them. Substitution of the 
ansatz 

f n (r,t) =5>„,-(i)/f (r^e-H 4 E ^ dt ' (14) 

3 

into Eq. ( |l0| ) yields the coupled differential equations 
^[a nj ff+a nj fP}e-*i:^'=o, (15) 
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where the dots above symbols denote time derivatives. 
Taking the positive scalar products of these equations 
with the state fg°\ utilizing the orthonormalization re- 
lations (|l3|), and solving for a n o, we find 

3 

In a similar manner, we derive the equations 
a nk = -J2a n je-*f> (t,) - Ek(t ' )]d Uf^ l/f )_ (16b) 

3 

for the coefficients corresponding to positive-norm states. 

In order to estimate the decay of a state jn \ we as- 
sume that at time t = its expansion coefficients are 
a n j(0) — S n j. Approximating the slowly varying scalar 
products and energy eigenvalues to be constant in time, 
and the decay to be negligible, such that we may also 
treat the a n j coefficients on the rhs of Eqs. ( |l6| ) as con- 
stants, we can integrate them to yield 



o,nk(t) — (e 



i)</W>±- 



(17) 



Above, the positive scalar product is chosen for k = 0, 
the negative otherwise, and we have denoted u> n k = 
[E n — Ef.)/h. Since the quasi-stationary states are or- 
thonormalized, the requirement of negligible decay thus 
implies 



1 



</fi/r> 



< i. 



(18) 



Essentially, this is the validity criterion of the adiabatic 
approximation for the TDHFBP equations of the dilute 
boson gas. 

Consider now the case of a vortex line precessing with 
frequency v w about a circular orbit of radius i? pr in a 
harmonically trapped condensate; for simplicity, we as- 
sume a trapping potential of the form V tI — ^maj^r 2 in 
cylindrical coordinates r = (r,0, z), and the vortex line 
to be directed along the z-axis. In view of the differences 
in the vortex-core structure between the BA and the self- 
consistent approximations, it is especially interesting to 
find out whether the lowest-energy quasiparticles, which 
constitute the major contribution to the noncondensate 
filling the vortex core, can follow the moving vortex line 
rigidly, i.e., adiabatically. In order to assess the validity 
of the criterion (|l^) for such states, we use the estimate 



</nv/r> : 



(19) 



where v is the velocity of the vortex line. This approxi- 
mation treats accurately the region in the vicinity of the 
vortex line, although it is exact only for a uniform vortex 
motion. Furthermore, supposing the precession orbit is 



not too near the condensate boundary, we may use the 
quasiparticle states of a system with a vortex located in 
the center of the trap to estimate the scalar products on 
the rhs of Eq. ([l9|) . Such a system is cylindrically sym- 
metric, and the quasiparticle eigenstates can be chosen 
to be of the form 



v g (r) 



(20a) 
(20b) 



where qe and q z are integer angular and axial momen- 
tum quantum numbers, respectively, and q denotes the 
complete set of quantum numbers for the states. Calcu- 
lation of the required matrix elements is straightforward 
for these states — the result is 



«V=v-</< J|V/«?>) 



dr 



dr 



±v! 



dr 



+(v'e ~ Qe)[(q'e + l)u*u g , ± (q' e - l)v*v q > 



(21) 



where v = |v| is the magnitude of the velocity of the 
vortex line, and the states are normalized according to 
J °° rdr ( | Ug | 2 ± |v 9 | 2 ) = 1, with the plus (minus) sign 
used for the zero-energy condensate state (other states). 
Equations (|l8|), (S9h, and (fjl|) finally yield the criterion 



v <C v, 



qq' 



lqq< 



(for all q) 



(22) 



for the velocity of the precessing vortex in order for the 
state to follow the vortex rigidly. 

We have numerically computed the adiabaticity ve- 
locities v qq i for the lowest excitations of a cylindrical 
condensate. The static HFB equations were solved self- 
consistently within the PA and its so-called Gl and G2 
variants []l2| , [l6|| , in order to find the quasiparticle am- 
plitudes u q (r), Vq(r), and the respective eigenenergies — 
for details of the methods used in the computations, see 
Ref. po| . In order to facilitate comparison with the 
vortex precession observations, we use parameter values 
which essentially correspond to the experiments reported 
in Ref. M. Especially, the radial trapping frequency 
was set to v r = uj r /2'K = 7.8 Hz, and the density of 
the trapped 87 Rb atoms was adjusted to yield a heal- 
ing length £ = (87rnoa) 1//2 ~ 0.7 /im at temperature 
T rj 0.8T B ec- Here no denotes the peak density of the 
condensate, and the condensation temperature T BEC « 30 
nK. 

In the experiments, the observed precession radii were 
of the order of i? pr ~ i?/3 ~ 10 /im, where R denotes 
the radius of the condensate JtJ . Bare-core vortices were 
observed to precess in the direction of the condensate 
flow with frequency v r * 1.8 Hz, which corresponds to a 
velocity v cxp = 2irR pr v r w 0.1 mm/s. This is to be com- 
pared with the computed velocities v qq i for the lowest 



3 



quasiparticle states with q z — 0, displayed in Fig. [I] [ pj| . 
Although v cxp < , we find the adiabaticity condition 
( p2| ) not to be fulfilled. This suggests that, due to the de- 
formation of the quasiparticle states, the noncondensate 
can not follow the vortex line rigidly at these velocities. 
Especially interesting is the smallest adiabaticity veloc- 
ity given by the decay of the so-called lowest core local- 
ized state (LCLS), which is the lowest excitation with 
(l8,qz) = (—1,0), and itself corresponds to the preces- 
sion of the vortex. The LCLS has a crucial role in the 
filling of the vortex core with noncondensate, which sta- 
bilizes the static vortex state [p|,[To|] . In fact, this state is 
almost solely responsible for the differences in the vortex 
structure between the BA and PA in the low-temperature 
limit. Deformation of the LCLS due to the vortex mo- 
tion thus implies crucial modifications for the vortex core 
structure. 

The given adiabaticity velocities also hold for the Gl 
and G2, since differences between the approximations 
turn out to be negligible in this respect. Computations 
with various parameter values also confirmed the valid- 
ity of the criterion (^2|) to be largely independent of the 
specific values of the trapping frequency, the density of 
the gas, or the effective interaction between the atoms. 
Essentially, the adiabaticity of the system is determined 
by the precession radius i? pr , via its proportionality to 
the velocity of the precessing vortex line. In addition, 
the adiabaticity could depend on the temperature: al- 
though we found the smallest v qq i to depend only weakly 
on temperature, the stationary PA predicts the preces- 
sion mode frequency and, hence, the precession velocity 
to have a strong temperature dependence p0| . 

In conclusion, we have derived a criterion for the va- 
lidity of the quasi-stationary approximation for a time- 
dependent mean-field formalism describing the dynamics 
of the condensate and thermal components of a dilute 
boson gas. Application of this criterion to a harmon- 
ically trapped Bosc-Einstcin condensate containing an 
off- axis, precessing vortex line is shown to yield for the 
vortex velocity a condition which is not fulfilled in the 
experiments conducted so far. Deformation of the vor- 
tex structure due to its motion is thus suggested to be 
at least partly responsible for the apparent discrepancies 
between the predictions of the stationary self-consistent 
approximations and the results of the vortex precession 
experiments. 
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